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Abstract 

Sutherland's theorem dictates that the contribution of the electromagnetic interaction to the 
decay process n — > 3vr° is neglected with respect to the one coming from the difference between 
the up and down quark masses. In the framework of chiral perturbation theory including virtual 
photons, we calculated the main diagram concerning the exchange of a virtual photon between 
two intermediate charged pions. The correction induced by this diagram on the slope parameter 
amounts to 17% of the correction induced by the pure strong interaction at one-loop level. If this 
result is maintained when considering all the diagrams at the chiral order we are working, we can 
say without any doubt that Sutherland's theorem is strongly violated. As a direct consequence, 
any determination of light quark masses from the present decay should take into account the 
electromagnetic interaction. 
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I. INTRODUCTION 



The decay rj — > 7r°7r 7r° is forbidden by isospin symmetry of the strong interaction. This 
symmetry is explicitly broken by electromagnetism on the one hand, and by the difference 
between up and down quark masses on the other hand. The decay amplitude takes then 
contributions from the electromagnetic fine structure constant af s and — m u . Suther- 
land's theorem states that the electromagnetic contribution is suppressed with respect to 
the other one 

HQ. 

As a consequence, the decay r] 3ir has been considered over the 
last thirty years as the golden process for the determination of the difference rru — m u by 
comparing the theoretical predictions with the experimental measurements (See [3| and ref- 
erences therein). However, the data for the slope in the corresponding Dalitz plot show a 
big discrepancy between the experimental value and the theoretical predictions based on 
chiral perturbation theory 4]. This fact motivated the scientific community to perform very 
accurate measurements at the experimental level ^-7]. On the other side, several attempts 
have been conducted at the theoretical level. For instance, a model-independent study rest- 
ing only on chiral symmetry and experimental data was performed in 8j. The authors of 
[2I used another approach based on data and analytic dispersive representation including 
first-order isospin breaking effects. Rescattering effects were considered in (9] through a 
modified non relativistic effective field theory beyond one loop including isospin-breaking 
corrections. It was claimed in the preceding citation that "the effect of photon exchange 
inside the charged pion loops on the Dalitz plot expansion is small, even on the scale of the 
other small isospin-breaking effects". We believe, and will partially prove, that this is not 
the case. Our statement rests on the observation that the re-scattering effects between two 
pions in the final state is dominated by the Coulomb interaction which, in two different (but 
related) processes, gives a relatively valuable correction 10Hl2| . Motivated by this observa- 
tion, we calculate in the present work the correction induced on the slope parameter of the 
decay 77 — > 3tt° by photon exchange between two intermediate charged pions. To do so, the 
corresponding Feynman integral has been reduced to a set of five Master Integrals following 
Tarasov's reduction algorithm [13I. [l4|. which is implemented in computer algebraic systems 



as a Mathematica package called Tarcer [151]. The Master Integrals are then obtained in 
dimensional regularization by solving hypergeometric differential equations as reviewed in 
. As a next step, we expand the decay amplitude around its value at the center of the 
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Dalitz plot and derive an analytic expression for the slope parameter. The value of the slope 
is then deduced in four dimensions by expanding around D = 4. We finally comment on 
our finding. 



II. THE DECAY AMPLITUDE 

The dynamics of the process 

V(p)^^( pi ) + n (p 2 ) + n (p 3 ) (1) 
is studied in terms of the Mandelstam variables 

s = (p~Pi) 2 , t=(p-p 2 ) 2 , u = (p-p 3 ) 2 } (2) 
subject to the constraint 

s + t + u = Ml + 3M^ = 3s . (3) 
The decay amplitude M. is defined as 

(vrVvr ! rf) = i(2rt) 4 6\ Pl + p 2 + p 3 - p)M , (4) 

and can be written by symmetry considerations like 

M(s,t,u) = M(s)+M(t)+M(u). (5) 

In the framework of chiral perturbation theory, we write the s-channel amplitude as 

Mis) = - B0{ ^ F T U) t 1 + + <W«) + firOO] + 8W) > (6) 

where S stT is of Q(p 2 ). S em of C(e 2 ), and S em of 0(e 2 p 2 ) in the chiral counting and have been 



calculated in jl7|, [l8|, and [8|, respectively. The correction 5 7 is of 0(e 2 p 2 ) and is supposed 
to be small compared to the others. The main diagram contributing to <5 7 is sketched in 



Fig. [TJ Working in the framework of chiral perturbation theory including photons 19( , we 
find that, up to polynomials in the kinematic invariants, the main diagram is given in terms 
of D-dimensional integrals of the form 

1 ff d D k 1 d D k 2 {k 2 ) i {k 2 y{q-k 1 ) k {q-k 2 y{k 1 -k 2 ) m 



7T 



D J J (k 2 - Ml){k\ - M%)[(ki - q) 2 - M%)[{k 2 - q) 2 - M%]{h - k 2 f ' 



(7) 
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Figure 1 : In this diagram q = p — p\ = p2 + P3 ■ 

Applying Tarasov's reduction algorithm, integrals ([7]) are reduced to a set of five Master 
Integrals, Ji, J 2 , JT, J 2 , and T 2 with 

1 f d D k 

J = J-f ft (9) 

J= J_ff fhfh (w) 

1 " ^ J J (k 2 - M*)[(ki - q) 2 - - h) 2 ' 1 ) 

\_ rr d D k x d D k 2 (k x ■ fc 2 ) 

We found that the contribution of the main diagram to S 1 takes the following form in D 
dimensions 



3F 2 M 2 (M 2 -M 2 o) S 3 (T 4 (s) 
(D _ 4KB -3) (4tt)^ (CiJi + C2 ' h + CsJT + + CsT2) 



(12) 



where 



and, 



4M 2 

4) = ll--- (13) 



ci = -(-384D 2 Mf + 1536DM* - 1152M 10 + 608D 2 sM^ - 2528DsM* 

+864sM* + 184£>VA^ - 624Ds 2 M^ + 2240s 2 M^ - 650L>VM 4 + 2810L>s 3 M 4 
-3408s 3 M 4 + 257L> VM 2 - 1299L>s 4 M 2 + 1564s 4 M 2 + 12L>V - 42Ds 5 + 48s 5 ) 

c 2 = 2(-576D 2 M* + 2880UM* - 3456M® + 1360D 2 sM^ - 6880.DsM£ 
+7776sM^ - 980L>VM 4 + 5168L>s 2 M 4 - 5568s 2 M 4 + 187L>VM 2 
-1213L>s 3 M 2 + 1356s 3 M 2 + 36£>V - 126L>s 4 + 144s 4 ) , 



c 3 = -sV(s)(56D 2 M° - 37QDM% + 416M° - QQD 2 sM 4 + 522DsM 4 - 600sM 4 

+D 2 s 2 Ml - 119Ds 2 Ml + 148s 2 M 2 + 12D 2 s 3 - A2Ds 3 + 48s 3 ) , 
c 4 = 3sV(s)M 2 (3s - AM 2 )(s - M 2 )(D - 3)(-ADM 2 + 8M 2 + 3Ds - 8s) , 
c 5 = -\92D 2 Ml + llh2DMl-lh3QMl + hmD 2 sMl-33l2DsMl 

+37\2sMl - 39QD 2 s 2 M 4 + 2568Ds 2 M 4 - 2640s 2 M 4 + 31D 2 s 3 M 2 
-489L>s 3 M 2 + 524s 3 M 2 + 36D 2 s 4 - 126L>s 4 + 144s 4 . 



III. THE MASTER INTEGRALS 

Three of the Master Integrals are given in terms of one-loop integrals T and J, 

T = -iM°- 2 T(l - D/2) , 
J = i jl -M° " 4 r(l -D/2) + a 

The remaining two integrals are found to satisfy the system of first order differential 
tions, 

d_h (5-2D)s + 2(3-2D)M 2 6(D - 2) 2(D - 2) 2 

ds s(AM 2 - s) 1 s(AM 2 -s) 2 s(AM 2 - s) ' 

(U2 8{D - 1)M 4 + 2{3D - 8)sM 2 + (D - 2)s 2 
ds ~ As(AM 2 -s) 1 

D-2 ( 8M 2 + s \ D-2 ( AM 2 + s \ ^ 
2s \AM 2 -s) 2 4s \AM 2 -s) 

Such a system is equivalent to a second order differential equation 

d 2 J, (D - 2f t2 3(4M 2 + Ds- 4s) dJ x 

ds 2 ~ 2s 2 (AM 2 -s) 2s(4M2-s) ds 

(D-3)[2(D-2)M 2 n + (D-A)s] 

2s 2 (AM 2 - s) 1 ' 

Differentiating once more, we obtain a third order differential equation 

d?Ji (D-A)(D-3) 28M 2 + 3Ds-18sd 2 J 1 
'ds' 3 ' ~ 2s 2 (AM 2 - s) 1 2s(4M2 - s) "ds 2 "" 
2D 2 Ml + sD 2 - 10DM 2 - 13s£> + 36s dJi 
+ 2s 2 (4M 2 - s) Is' 

Making the change of variable, 



In 



a 



VK 



we find that J\ satisfies the hypergeometric differential equation 



x 2 (l +x) 



d 3 J x 
dx 3 



+ 



3(D-6) 7 

— x x 

2 2 

(D-A)(D-9) D(D — 5) 
-x — 



2 4 
The general solution of such an equation is 

(3-D), (4-£>)/2,l 



dx 2 

dJi (3-D)(A-D) 



+ 



Ji = 0. 



(22) 



<A — 3-^2 



D/2, (5 - D)/2 



— x 



A 



+ (-x)-^l\F 2 



(8-3D)/2,(3-D),(A-D)/2 



; -x 



B 



(4-D)/2,(7-2D)/2 

(2D - 3)/2, (£> - l)/2 ' 
In order to determine the integration constants A, B, and C, we notice that 



+ (-x)^- 3 )/ 2 3 F 2 



C. 



(23) 



■Mo) = 



D — 2 



rril 
■*■ 1 



(24) 



2M2 VD-3, 

is regular for 2 < D < 4. It follows that B = C = since the terms (— x)~( D ~ 2 ^ 2 and 
(— x)( D ~ 3 ^ 2 are divergent for x = and in the range in question. Moreover, 3 F 2 is equal to 
one for x = 0. This means that A = Ji(0) and 



J 1 = ^,(DZ*M 



2M 2 \D-3 
Having J\, we obtain J 2 from 



(3-D), (A-D)/2,l 
D/2,(5-D)/2 



-x 



rp2 



(25) 



1 2 2Ml 
h ~ 3 + 3(D^2) 



x(l + x)^- M ^ 2D -^t 3 " 2D] ^, 
V ' dx 3(D-2) 



(26) 



by the use of 



dx 



3-^2 



(3-D), (4-£>)/2,l 
D/2,(h-D)/2 



; -x 



2(D-3)(D-4) 
D(D-5) 



■ 3-^2 



(4-£>),(6-D)/2,2 
(£) + 2)/2,(7-D)/2 



; -x 



(27) 



We get 



J. 



2(2£> - 5)x + 3 - 2D 



6(D - 3) 



■ 3 r 2 



(3-D),(A-D)/2,l 
D/2, (5 - D)/2 



; -x 
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2(D-4) , 



(4-C),(6-D)/2,2 
p + 2)/2,(7- J D)/2' 



-i^T 2 . (28) 



IV. THE SLOPE PARAMETER 

The amplitude is expanded around its value at the center of the Dalitz plot as 

\M(s,t,u)\ 2 = \M(s ,s ,s )\ 2 {l + 2az} , (29) 

where 



z = 



2M 2 (M V - 3M w o) 2 
We write the slope parameter as 



{(s - s ) 2 + (t- s ) 2 + (u- s ) 2 } . (30) 



M 2 ± - M\ . . 

a = a stv + a cm + a 7 + -—f -^a e m • (31) 

The contribution of the main diagram to the slope parameter is 

M 2 

a 7 = -^(M v - 3M,o) 2 Be 6" (s ) . (32) 

y 

We can express a 7 in terms of the five Master Integrals thanks to the differential equations 
satisfied by the latter. Up to now, we have considered only the differential equations satisfied 
by J i and J 2 . It is easy to check that the one- loop integral satisfies the following differential 
equation 

dJ _ (D-A)s + AM 2 D-2 

ds 2s(s-AM 2 ) s(s-AM 2 ) ' 1 ' 

Differentiating 5 1 twice with respect to s and using the preceding relation at each step, we 
obtain 



a is M 2 (M v -3M n o) 2 1 
W8F 2 M 2 (M 2 -M 2 ) sla%s ) * 

^Re(d 1 J 1 + d 2 J 2 + d 3 JT + d 4 J 2 + d 5 T 2 ) s=s , (34) 



(D - A)(D - 3) (Ait) - 1 v 11 z z ° " Js =*o 

with 

di = 6144(D 4 - 16D 3 + 86D 2 - 176D + 105)M^ 4 - 512(25£> 4 - 360Z) 3 + 1511D 2 
-2214D + Q21)sM^ + 128(48L> 4 - 374L> 3 + 257L> 2 + 3240L> - 5856)s 2 M 40 
-32(50L> 4 + 2009£> 3 - 17080D 2 + 43918D - 36964)s 3 M| + 8(397L> 4 - 684£> 3 



-15994D 2 + 65826D - 69384)s 4 M^ + 2(201£> 4 + 6971D 3 - 60836D 2 
+ 161172D - 139248)s 5 M 4 + (-1697D 4 + 16825D 3 - 62006D 2 
+101120D - 61664)s 6 M 2 + 6D(-2D 3 + 11D 2 - 22D + 16)s 7 , 

d 2 = 2[-9216(L> 4 - 17D 3 + 101D 2 - 247L> + 210)M 1 . 2 + 256(103L> 4 - 1527£> 3 
+7880D 2 - 17199D + 13446)sM 10 - 64(454L> 4 - 5634£> 3 + 25631D 2 
-50730D + 36744)s 2 M^ + 48(398L> 4 - 3551D 3 + 12214£> 2 - 19708D 
+12600)s 3 M^ + 4(-2491£> 4 + 18774D 3 - 43814D 2 + 26916D + 11088)s 4 M 4 
+ (2971£> 4 - 30543L> 3 + 114482D 2 - 188208D + 115296)s 5 M 2 
+18D(2D 3 - 11D 2 + 22D - 16)s 6 ] , 

d 3 = s 2 a 2 (s)[38A(10D 3 -77D 2 + 207D - 172)M^° + 32(18 J D 4 - 320D 3 
+1733D 2 - 3751D + 2764)sM^ - 8(187L> 4 - 1741L> 3 + 5696D 2 
-7898D + 3816)s 2 M^ + 6(209£> 4 - 1951D 3 + 5866D 2 - 6400L> 
+1568)s 3 M 4 + (-325D 4 + 3953D 3 - 16398D 2 + 29008D 
-18848)s 4 M 2 + QD(-2D 3 + 11D 2 - 22D + 16)s 5 ] , 

d 4 = -l2M\s 2 a 2 (s)[5l2(D 2 -5D + Q)Ml° + 128(D 3 -IID 2 

+36D - 36)sM^ + 16(D 4 - 16D 3 + 91L> 2 - 216D + 180)s 2 M° 
-8(5L> 4 - 38£> 3 + 92L> 2 - Q5D - 12)s 3 M 4 + (D — 3) 2 (33D 2 
-Q8D - 64)s 4 M 2 - 3(D - 3) 2 (3D 2 - UD + 16)s 5 ] , 

d 5 = -3072(L> 4 - 18L> 3 + 115D 2 - 306L> + 280)M 1 . 2 + 256(37L> 4 - 615£> 3 + 3401D 2 

-7860D + 6444)sM 10 - 64(183L> 4 - 2517D 3 + 12185D 2 - 25152D + 18772)s 2 M® 
+16(568L> 4 - 5294D 3 + 18567D 2 - 30394D + 19784)s 3 M^ + 4(-1386£> 4 + 10909D 3 
-25796D 2 + 15140D + 7808)s 4 M 4 + 5(359£> 4 - 3883D 3 + 14986D 2 - 25136D 
+ 15648)s 5 M 2 + 18D(2D 3 - 11D 2 + 22D - 16)s 6 . (35) 

V. RESULTS AND CONCLUSIONS 

We first expand 5 1 around D = 4 and then use the following numerical values 

"fe = 7^77777 > ( F -> M -° = M -' M v) = ( 92 - 42 - 139 - 57 > 547 - 3 °) MeV • (36) 
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We found that the amplitude, <5 7 , possesses a pole at the lower bound of the allowed kine- 
matical region, that is, for s = AM 2 . At the upper bound, s = (M v — M n o) 2 , the calculated 
electromagnetic correction amounts to 3-4% of the Tree level amplitude. This points out 
that the present correction can not be simply neglected as dictated by Sutherland's theorem 
and as is widely believed by physicists. 

We next do the same for the slope parameter and obtain 



This value is to be compared with the one-loop strong and electromagnetic corrections, 



respectively. The first conclusion we draw comes from the comparison between the electro- 
magnetic corrections at one- and two-loop levels. Although the two-loop calculation does 
not include the contribution of all the diagrams, it is three times bigger than the one-loop 
calculation. Since the one-loop correction is originating from the mass difference between 
charged and neutral pions, and the two-loop correction considered here is due to photon 
exchange between charged intermediate pions, we can claim that the effect of soft virtual 
photons is more important than the one for hard virtual photons. Obviously, this conclu- 
sion does no more hold if some cancellations occur between the various contributions if we 
consider all diagrams with a photon exchange. This point should constitute a sufficient 
motivation for calculating the full amplitude including the exchange of one photon between 
charged intermediate light mesons. 

We draw the second conclusion by comparing the size of the electromagnetic correction 
to the size of the one-loop pure strong correction. If we consider only the contribution 
of the photon exchange diagram, we find that the former amounts to 17% of the latter. 
Adding both one- and two-loop corrections, we find that the electromagnetic interaction 
contributes about 11% of the pure strong interaction at one loop. If we add both strong 
and electromagnetic corrections, at the one-loop order for the first, and the two- loop order 
for the second, we obtain for the slope parameter 



= 0.0029 . 



(37) 



a str = 0.0179, 



a em = -0.0011 , 



(38) 



a = 0.0179 + -0.0011 + 0.0029 




+0.0198, 



(39) 



0(p 4 ) 0(e 2 p 2 ) 0(e 2 p 4 ) 
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to be compared with the most recent experimental value [7J, 

a exp = -0.0301 ± 0.0035i°;°o35 , (40) 

noting this endless sign discrepancy between the chiral perturbation theory prediction and 
observation. 
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